Angular conductance resonances of quantum dots non-collinearly coupled to 

ferromagnetic leads 
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The zero bias conductance of quantum dots coupled to ferromagnetic leads is investigated. In 
the strong coupling regime, it is found that the conductance is a non-monotonic function of the 
angle between the magnetisation directions in the two contacts. This behaviour is an effect of the 
presence of the leads which induces an angle dependent spin split of the quantum dot states, and 
spin flip transitions between the quantum dot states whenever the magnetisation directions of the 
leads are non-collinear which enhances the current density at the chemical potential. In the weak 
coupling regime, the system reverts to normal spin valve character. 
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PACS numbers: 72.25.Mk, 73.63.Kv, 73.23.Hk 



Due to the fundamental complexities and far-reaching 
technological possibilities, the interest in spin-dependent 
transport in mesoscopic systems remain as high as ever. 
Since the first measurements of giant magneto resistance 
(GMR) in Fe/Cr magnetic super-lattices 0, magneto- 
transport has been studied for normal or ferromagnetic 
metallic islands 2], and spin-dependent transport from 
ferromagnets through quantum dots (QDs) [3, Q and 
molecules 0. Recently, Kondo physics of QDs weakly 
coupled to ferromagnetic leads have been extensively 
studied 0, 0, H|, showing that the suppression of the 
Kondo resonance due to the spin-polarisation of the leads 
may be restored by application of an external magnetic 
field. 

In this Letter, effects of strong Coulomb interactions 
on the linear conductance is studied, in a single level 
QD coupled to ferromagnetic leads in the absence of 
external magnetic fields. The conductance of the sys- 
tem is predicted to display a huge complexity, for non- 
collinear alignment of the magnetic contacts. Due to the 
strong on-site correlations, the localised states undergo 
a spin-split depending on the non-collinear magnetisa- 
tion of the leads, which in combination with effects from 
spin flip transitions, result in a non-monotonic conduc- 
tance in the strong coupling regime. The QD system is 
analysed in terms of non-equilibrium many-body opera- 
tor Green functions (GFs) [!J, 113 which is motivated for 
three reasons, namely, 1) the on-site Coulomb repulsion 
is the largest energy scale of the system, 2) a renormali- 
sation of the QD level similar to the the scaling relation 
found in |(J is included in the QD GF 0, which cannot 
be obtained within traditional standard methods, and 3) 
the theory is valid in the whole range from the weak to 
the strong coupling regime ^l|. In the present case the 
Kondo effect is suppressed by the ferromagnetism in the 
leads 0, Q, Q , and can therefore be neglected. 

Consider a single level QD with a bare quantum level 
which is spin-degenerate, in the atomic limit, at the 




FIG. 1: Quantum dot in non-collinear coupling to the fer- 
romagnetic leads. The co-ordinate system (left) defines the 
global reference frame. The magnetisation directions of the 
leads (arrows) enclose the angle cf>. 



energy e , and that the on-site Coulomb repulsion is 
given by U. Hence, the energy of the QD is given by 
Hqd = J2a £ od^d a + Un^ni, where d\ (d a ) creates (an- 
nihilates) an electron in the QD, whereas n a = d\d a , 
and a =1", \ is the spin projection in the global refer- 
ence frame. Since the Coulomb repulsion is the largest 
energy scale, the QD is rewritten in terms of its eigen- 
states, e.g. \p), p — 0,(7,2, where |2) = | fj,), hence 
KQD = E p E P hP, where W = X™, {X™ = \p){q\), and 
Eq = 0, E a = Sq, and E 2 — 2e + U are the energies of 
the empty (|0)), singly (\a)) and doubly (|2)) occupied 
states, respectively. 

The leads are modelled by H L / R = Y^ka&L/R^Aa^ 
where ct (k), creates (annihilates) an electron in the 
left /right (L/R) contact at the energy e^. The mag- 
netisation in the left contact coincides with the z-axis of 
the global reference frame, see Fig. ^ hence the tun- 
nelling interaction between the QD and left contact is 
given by Htl = Y.kveL,a( v k<?{d (y ) a cl <y X a + H.c), where 
£ a (d CT ) a X a = (0\d ry \a)X Oa + {a\d a \2)X s2 , and a is the 
opposite spin of ct. In the following, it will assumed that 
the Coulomb repulsion U is sufficiently large so that the 
doubly occupied state can be neglected, for briefness. 
Thus, EaK) a ^ a = X 0a , since (0\d a \a) = 1, which 
yields Htl = E toei (^4^ 0<r 



H.c). However, the 
doubly occupied state is easily included by a straight for- 
ward generalisation of the results presented below, and 
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has been included into all numerical examples given. The 
magnetisation in the right contact is rotated by the angle 
in the global xz-plane, see Fig. Therefore, the tun- 
nelling between the QD and the right contact is modelled 
by 
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where the vector 
matrix 



(4 



fe+ 



and the spin rotation 



cos 0/2 sin 0/2 
— sin 0/2 cos 0/2 



Here, the spin indices s = ± are used because of the ro- 
tated magnetisation in the right lead. Letting the mag- 
netisation in the left lead coincide with the global z-axis 
is not a restriction since the magnetic properties of the 
QD depends on the angle between the magnetisation di- 
rections in the leads. 

The total Hamiltonian for the system is given by 
H = ^2 a=LtR CH a + Hto) + Ti-QD (stray fields from 
the leads are neglected). The spin-dependence of the 
leads is modelled by the hybridisation functions T"(ui) — 
Kl 2 %-^) = 2ttK| 2 ^H, where t& = w« 
and is the spin-dependent density of states in the leads 
which is assumed to be constant. The spin-dependence of 
r£ is parametrised in terms oip a = (T* -T*)/ (T* + 
letting = r (l ± p a ), where 2r = T" + T". By 
thisprocedure no essential physics is lost, as discussed 
in jjj. In terms of the spin-dependent parameters p a , 
the coupling matrices to the left/right leads become 
r L = r diag{l + pl, 1 — Pl} (diagonal matrix) and 



t r _ r / 1 + PR cos p R sin 
1 p r sin 1 — pr cos ( 



A straight forward derivation shows that the zero bias 
conductance for this system can be written as 

G(p,cf>) = ^ / tr [r L G r (co,p^)T R G a (co,p^)] 

x^cosir 2 (/?^^ (i) 

for low temperatures (/3 _1 = /cbT), where p is the 
chemical potential of the system. In this expres- 
sion, G r / a (w,/i,0) = D r / a (w,^i,0)P(/i,0) is the re- 
tarded/advanced Green function (GF) of the localised 
states in the QD, where P(p, 0) and D r / a (o;, p, 0) de- 
note the spectral weight and locator, respectively. The 
locator carries the local on-site properties of the GF, that 
is, the positions and widths of the poles. More details on 
the definition of the many-body operator GFs and the 
diagrammatic technique employed here can be found in 
0- fiol El ■ Following the steps in [|J and including the 



dependence of the coupling to the right, the QD GF is 
derived to (details will be published in [lj]) 

[ul - A(/i, 0) - £*-/«(//, 0)]C-/<>, p, 0) = P(p, 0), (2) 

where / is the identity matrix, whereas the transition 
energy matrix A(/i, 0) is renormalised by kinematic in- 
teractions between particles in the different localised 
states due to the presence of the de-localised electrons 
in the contacts. Here, the locator is identified by 
D r /«(u;,/i,0) = [luI - A(pi,0) - E^Cw.Ai,^]" 1 . The 
transition energy matrix is given by 



A(M, 



A°+ £ 



a=L,R 



f a (e) - f(w) 



x^r Q [-2ImD>, M ,0)K^^. (3) 

Z7T ZTT 

where A° = (eo — Eq)I is the bare energy matrix, 
f a {uj) — f(tu — p a ) is the Fermi function (pl/r is the 
chemical potential of the left/right lead; here p a = p), 
and a x is the cc-component of the Pauli spin vector. 
The self-energy in Eq. (0) is given by X r / a (/i, </>) = 
=FiP(^, 0) J2 a T Q /2, since the real part is negligible due 
to the large conduction electron band width in the leads. 
Apart from the renormalisation of the transition ener- 
gies, Eq. @ reduces to the result by Varma and Yafet 
|l5| in the non-magnetic case for large U . 

It should be emphasised that the QD GF has to be 
self-consistently calculated, for each point in the param- 
eter space (p,(j),T, A ,T ,pl,pr), along with Eq. © 
and evaluation of the end-factor P. The latter quan- 
tity is defined by P aa , = (T{A 0cr , X 17 ' }) = 5 a<T ,N Q + 
N a i a , that is, as a sum of the population numbers of 
the states involved. These are calculated by No — 
-Im£ a /G>(w)dw/(27r), N a = hn J G„{w)du /(?*), 
and Ne'e = —i J G^ a ,(u>)dLu/ (2ir) (note the reversed or- 
der of the spin indices), where G < / > = G r V < / > G°, 
Gc7^ y and G^T are the diagonal and off-diagonal com- 
ponents of G</>, and V< = i[f L T L + f R T R ] and 
V> = V< - i[T L + r R ], The self-consistent calculations 
are subject to the condition Nq+J^c = 1, which satis- 
fies the requirement that the integrated total QD density 
of states is unity, e.g. Im f tr [G < — G > }du)/(2ir) = 1. 
In this fashion effects from the spin-flip transitions are 
taken into account by self-consistently solving the full 
2x2 matrix equation, e.g. Eq. JQ), giving non- vanishing 
off-diagonal components whenever < 4>/tt < 1. 

The renormalisation of the transition energies, Eq. Ip^jl. 
arise due to kinematic interactions between particles in 
the localised states induced by the presence of the de- 
localised electrons in the contacts. The magnetic prop- 
erties of the contacts are transferred into the QD by this 
renormalisation, since the energy for the spin a state are 
influenced by the properties in the spin a channel of the 
system. This is easiest seen for collinear leads, since then 



FIG. 2: (Colour online). Renormalised transition energies 
as function of tf>. Here p L/R = 0.5 and {A , U, W, T}/To = 
{0,50,2000,0.17}. 



E ct5 , P aS = 0. Assuming zero width of the locator in Eq. 
(J3J) and large band-width 2W 3> \p — A a \ of the conduc- 
tion bands, result in [9j 




A„ = A 



a=L/R 



A, 



W 



(4) 



This expression is consistent with the results in 0, 0] 
for the non-magnetic case, as well as with the scaling 
equation reported in 6] for the magnetic case. The log- 
arithm is negative which leads to A a < A . By the 
replacement A a — > A on the right hand side of Eq. 
(QJ, e.g. the first iteration in the self-consistent calcula- 
tions, one finds the difference A| — Aj_ = — (T /n)(pL + 
Pncoscf)) log \{n — A°)/W\, in agreement with Con- 
sequently, this difference is positive (negative) for p a < 
(p a > 0), i.e. for parallel (0 = 0) magnetic leads 
such that > T^, (Fj* < T"), while for it vanishes 
for anti-parallel (<f> = 7r) leads such that pl = PR- The 
renormalised transition energies in these two limits can 
be viewed in Fig. EK^/tt = 0, 1). The result in Eq. © 
provides a continuous connection for the transition ener- 
gies between the two collinear configurations of the leads, 
where the maximal and minimal spin split in the QD are 
given for the parallel and anti-parallel cases, respectively, 
as shown in Fig. El The characteristics of the transition 
energies A a vary only slowly with p, and therefore the 
plots in Fig. [3 will be used independently of p. 

For non-collinear leads, e.g. < <fi/ir < 1, the physics 
of the QD becomes more subtle, however, since the off- 
diagonal elements of T R are finite, which leads to that 
the spin- flip transitions acquire non- vanishing weights, 
e.g. N a9 ^ 0, hence, the off-diagonal components of 
the QD GF G aa <(ioj) 7^ 0. Then, it is expected that 
both spin projections of the local density of states (DOS) 
p a = Im[G<— G>]/(27r) = — ImG a /ir become mixtures of 
one another in the sense that both p-\ and pi are peaked 
around both A-f and Aj, since for instance 



{uj~A- a - s;)p g + (A a9 + S^ g )P ge 
(w - zl){uj - z r 2 ) 



FIG. 3: (Colour online), a) Contour plot of the conductance, 
< G(p,0)/G o < 1/2, where Go = 2e 2 /h, as function of 
the chemical potential, p, and angle, <j>. b) The conductance, 
G(p, 4>)/Go, for various values of the chemical potential indi- 
cated in panel a) . c) Current density for different angles <f> at 
/i/To = —2.06 (straight solid). The second panel from above 
also shows the spin projected local DOS. Here Pl/r = 0.85 
and A /To = (straight dashed). 



where the complex roots are given by 



1,2) 



zl = lim(tr D r,_1 + (-l)Vtr 2 D r .- 1 - 4 det T>^)/2, 

and A gS , E ct5 . ^ 0. An example of the local DOS in this 
case is plotted in Fig. c) for <j)/-K — 2/3, where it is 
readily seen that both p^ (dotted) and pi (dashed) are 
doubly peaked, having finite densities in between. 

In the remainder of this Letter, I discuss the effects 
of the (j) dependent renormalised transition energies and 
QD GF on the zero bias conductance. The contour 
plot in Fig. a) displays G(p, <j>) as the systems is 
swept from weak coupling (p — max{A|,A|} ^ To) 
to strong coupling and then back to weak coupling 
(p — min{A|, A|} <C —To), explicitly showing the com- 
plexity of the expected conductance in different regimes. 
Note that the character of the conductance peaks resem- 
bles the shape of the transition energies in Fig. [21 In 
the weakly coupled regime, the QD effectively behaves as 
a spin-valve, where the conductance G(p, <p) is a mono- 
tonic function of 0, see Fig. 01b) (dotted, dash-dotted), 
in agreement with jj| . This is expected since the minimal 
distance min CT \p — A a \ grows with 0, which then leads to 
a reduction of the electron density in the QD around p. 

In contrast, when p lies in the vicinity of one or 
both transition energies (strong coupling), the quali- 
tative character of the conductance is dramatically al- 
tered. In the parallel configuration (0/7T = 0) where 
the spin split of the QD transition energies is maximal, 
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FIG. 4: (Colour online), a) Average spin-polarisation, Nf — 
Ni, and b) occupation, Nf + N±, of the QD as functions of 
the chemical potential and angle. System parameters as in 
Fig. [3] In b), the dots indicate the minimum occupation at 
the corresponding value of the chemical potential for rotation 
angles < 4>/n < 1. 

the two spin projected states, peaked at Af ^= Af, are 
separated roughly by To(pl +pn)/Tr, and with widths 
P a T [2 ± (p L +p R )]/2, since E ctS , P aS = 0. For finite 
</>/-7r < 1, the two spin projected states become inter- 
mixed due to spin flip transitions, i.e. T> aSl P aS ^ 0. 
This provides an enhancement of the electron density 
around /i which then is reflected in the current density 
j(uJ,ti,(p) = tr T L G r {uj,ii,<f>)T R G a {uj,ii,4>), see Fig. 
c), which displays j(u>,[j,, <fi) for various rotation angles. 
As seen in this plot, the amplitude of the current den- 
sity at /J, /To increases for growing rotation angles up to 
a threshold angle (f>*/Tr, < <fi* /tt < 1, and thereafter 
decreases. 

The effects of the </> dependent renormalised transition 
energies on the QD spin-polarisation, Nf — Nf , and occu- 
pation (valence) , Nf + Nf , provides as additional under- 
standing of the conductance. The plots in Fig. ^display 
a) the calculated spin-polarisation and b) the occupation 
of the QD in the (/i, 0)-plane, and it is clearly seen that 
the QD is strongly spin-polarised for all values of the 
chemical potential between the QD state energies, i.e. 
mhvlACTo} < M < max {A CT o}, f° r = 0. The high 
QD occupation and large spin-polarisation in this regime 
reveals that there is only a small fraction of the local DOS 
available for transfer of electrons through the QD, since 
the spin-polarisation of the QD is opposite to that of the 
leads. As expected from the discussion of the spin-split 
of the QD state energies, the spin-polarisation decreases 
and eventually vanishes as 4>/-k — > 1. Due to the spin-split 
of the QD states, it is also clear that the spin-polarisation 
in the regime /x — mino-jAo-o} < is smaller than in the 
regime /i — maxo-jA^o} > since the QD is almost empty 
in the former regime whereas it is almost fully occupied 
in the latter, c.f. Fig. Elb). As is signified by the marks 
in Fig. 21b), the occupation of the QD clearly has a non- 
monotonic 4> dependence in the strongly coupled regime 
(mincrlAcro} < M < m ax CT {A CT0 })- Although the char- 
acteristics of the occupation and the transport proper- 
ties not are completely correlated throughout the whole 
(p, 0)-plane, a closer analysis of the QD occupation, in 



the strongly coupled regime, reveals that its minimum 
occurs at the position of maximal conductance. 

In conclusion, it is predicted that the linear conduc- 
tance of a QD strongly coupled to ferromagnetic leads is 
a non-monotonic function of the angle between the mag- 
netisation directions of the leads. The predicted charac- 
ter is an effect of an angle dependent spin split of the QD 
state, induced by the presence of the leads, and spin flip 
transitions between the QD states whenever the mag- 
netisation directions of the leads are non-collinear. In 
accordance with previous considerations Q, the system 
reverts to normal spin valve character in the weak cou- 
pling regime. 
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